Spin analog of the controlled Josephson charge current 
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We propose a controlled Josephson spin current across the junction of two non-centrosymmetric 
superconductors like CePtaSi. The Josephson spin current arises due to direction dependent tun- 
neling matrix element and different momentum dependent phases of the triplet components of the 
gap function. Its modulation with the angle £ between the noncentrosymmetric axes of two super- 
conductors is proportional to sinf;. This particular dependence on £ may find application of the 
proposed set-up in making a Josephson spin switch. 

PACS numbers: 74.50.+r, 72.25.-b, 74.20.Rp, 74.70. Tx 



Traditionally Josephson junctions [1, 2] in supercon- 
ductors draw interest both scientifically and its applica- 
bility in making devices. With no exception, they have 
also been studied in unconventional superconductors like 
spin-singlet cuprate [3] and spin-triplet Sr 2 Ru04 [4] su- 
perconductors. However no Josephson junction between 
nonmagnetic superconductors is known to generate spin- 
polarized current. The purpose of this letter is to the- 
oretically show that the direction dependent tunneling 
matrix element across the junction between two recently 
discovered [5] non-centrosymmetric superconductors like 
CePtaSi, leads to tunneling of both spin-singlet and spin- 
triplet Cooper pairs. As a consequence, nonvanishing 
spin- Josephson current is viable along with the usual 
charge- Josephson current. This novel spin- Josephson 
current depends on the relative angle £ between the axes 
of non-centrosymmetry and fiR in the left and right 
side of the junction respectively. This angular depen- 
dence may be used to make a Josephson spin switch. 

The normal state Hamiltonian [6, 7] for the electrons 
in a band of a lattice without inversion symmetry is 
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where electrons with momentum k and spin s (=t or |) 
are created (annihilated) by the operators c ks (cks), and 
£k is the band energy measured from the Fermi energy 
tp. The second term in the Hamiltonian (1) breaks par- 
ity symmetry as g_k = — <7k for a non-centrosymmetric 
system. For a system like Heavy fermion compound 
CePtaSi which has layered structure, Ho is considered 
to be two-dimensional. For such a system of electrons 

i 2 

with band mass m, £k = ^ — e F and — arj^ where 
r/k = ri x k, i.e., the spin-orbit interaction is of Rashba 
type [8] and a is the Rashba parameter. Here h repre- 
sents the axis of non-centrosymmetry which is perpen- 
dicular to the plane of the system. Due to the breaking 
down of the parity, spin degeneracy of the band is lifted; 
by diagonalizing Ho, one finds two spin-split bands with 
energies £kA = £fc + Aa|k| where A = ± describes helic- 
ity of the spin-split bands. Therefore in the diagonalizcd 
basis H (1) becomes H = J2k.x=± £kAC kA Ck A , where 



c~kA = (ck| — *A exp (i0k)cki) / \/2 is the electron destruc- 
tion operator and c tA = ^c k | + iXcxp (— i<^k)c k jj /V2 is 
the electron creation operator in band A with momen- 
tum k whose orientation with cc-axis is 0k- The density 
of electronic states at Fermi energy in these bands may 
be found as ^ = ^ (l — Xma/y^kp + m 2 a 2 ^j , where 

kp = \J2mep is the Fermi momentum. 

Band structure calculation [9] on CePtaSi reveals that 
the energy difference between two spin-split bands near 
kp is 50-200 meV which is much larger than the super- 
conducting critical temperature, ksT c s» 0.06 meV [5]. 
The formation of Cooper pairing between electrons in dif- 
ferent spin-split bands may thus be ignored. The Hamil- 
tonian for these superconductors may then be written 
as 
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where Cooper pairs are only between intraband electrons. 
Therefore the normal and anomalous Green's func- 
tions are obtained respectively as Q\{k.,ie n ) = —(ie n + 
&x)/(el + & + l A kA| 2 ) and T x (k,ie n ) = A kx /(e 2 n + 
£ kA + |AkA| 2 ), where e„ is the fermionic Matsubara fre- 
quency [13]. The superconducting order parameter AkA 
obeys the symmetry [9, 10]: A_kA = —AkA- We consider 
AkA = Afe^Ak, i.e., the angular dependence of k on the 
order parameters of two bands are assumed to be same. 
Apart from the overall phase rigidity angle of the su- 
perconductor, there may be a relative phase difference 
between superconductors of two bands: Afe + = |Afe + |e 
and A fe _ = |A fe _ \e l(e+e \ Reverting Hi (2) to spin-up 
(-down) basis, we find 
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A k , u = 

iexp[2i^k](A fe+ + Afe_)A k , and A k , u = -A k ,xT = 
| exp [i(/>k](Afc + — Afc_)Ak are the different compo- 
nents of pairing potential Ak, S s' between electrons with 
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FIG. 1: Schematic picture of two non-centrosymmetric su- 
perconductors, denoted by L and R, coupled at i = xq. 
Their orientations are characterized by the axes of non- 
centrosymmetry, til and ha. 



spins s and s' , and Tk.R = ia\k\ exp [— z</>k] is the 
Rashba spin-orbit coupling potential. The Hamiltonian 
(3) has been expressed in the basis such that \I/ k = 
(<vp c -kT> c kj.i c -ki)- Since A kjU = -A kjiT , there is no 
triplet component with zero projection along the spin- 
quantization direction; pairing between electrons with 
unequal spins entirely gives rise to singlet component. If 
we choose Ak = — zexp [— ifa], the triplet component of 
the pairing may be expressed as A^-.t = (dk ' <r)icr y with 
<ik = jjkT + Afe_)rjk. In other words, the only sta- 
ble [7] spin triplet component whose dk || corresponds 
to this Ak- Therefore Ak.ff and Ak,j.x will have equal 
but opposite momentum-dependent phase. Accordingly 
the singlet component of the pairing potential becomes 
A k ,s = 5(A fc+ - A fe _). If |Afe+| ^ |A fe _|, the admix- 
ture [6, 7] of singlet and triplet pairing takes pace. Re- 
cent observation [11] of Josephson current in the junc- 
tion of CePtsSi and s-wave superconductor suggests the 
existence of spin-singlet order parameter, while larger 
[5] upper critical field H c i seems to suggest that the 
spin-triplet pairing occurs. One normally assumes the 
superconducting order parameter to be independent of 
the magnitude of momentum, i.e., is fc-independent. 
The anomalous Green's function J r ss /(k., ie n ) in the up- 
down basis is related to T\ as T-\\ — + J 7 -), J~a = 
\e 2 ^{T+ + TJ), and T n = -T^ = - TJ). 
Apart from the applicable momentum-dependent phases, 
triplet (singlet) component of T ss i is the addition (sub- 
traction) of anomalous Green's function of two spin-split 
bands. 

We now consider Josephson tunneling between two 
such non-centrosymmetric superconductors as depicted 
in Fig. 1. The Hamiltonian for the system then reads 
H = Hl + Hr + Ht, where Hl and Hr are the 
bulk Hamiltonian of left and right side of the junction 



respectively and Ht describes the tunneling between 
these two sides. The Hamiltonian Hl is described by 
H\ (2) and Hr is also defined equivalently with the 
change in the notation of momentum k — * p to dis- 
tinguish each side of the junction. Further, each pa- 
rameter in the left (right) is denoted by superscript 
or subscript L (R). The tunneling Hamiltonian reads 

H t = Ek P ,s (^kpCk.Cps + T^pC+sCks) , where T kp is the 
tunneling matrix element for an electron with momentum 
p to tunnel from right side to left side with momentum k. 
Time-reversal symmetry of Ht suggests T_k,- P = T£ p . 

In a charge-tunneling process, number of electrons 
in either side of the junction changes with time, e.g., 
the rate of change in the number of electrons of the 
right side of the junction Nr = ^ p s c ps c ps reads as 

Nr — i[Ht, Nr]. Therefore the tunneling charge- 
current becomes I c (t) = — e^iV^(t)^: one part corre- 
sponds to quasi-particle tunneling and the other part 
is due to process of Josephson tunneling [12, 13], i.e., 
tunneling of the Cooper pairs. We are herewith in- 
terested in the Josephson charge current defined as 
I°j(t) = 2e3 [exp (-2ieVt) <f> c ret (eV)} at time t and for 
an applied bias voltage V across the junction, where 
$ c (iuj m ) = - f£ dr e luJmT (T t A(t)A(0)) with imaginary 

time r, A(r) = J2kp S T *p c 1.s( t ) c ps( t ')' tne bosonic Mat- 
subara frequency [13] uj mi and inverse temperature /3. 
Instead of the axes of noncentrosymmetry of two sides of 
the junction being parallel, general axes Hl and ur will 
lead to r/k = Ul x k and i) p = fiR x p. In that case, 
Josephson tunneling occurs in both singlet and triplet 
channels. Separating these channels, we find 
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where 
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XT = ^2 AkAp-F^k, ze„)J r A /(p, ie n - iuj m ) (6) 
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represent singlet and triplet contributions respectively. 
While the direction independent electronic tunneling 
probability |7k p | 2 = |T| 2 leads to tunneling of sin- 
glet component only, the triplet components of Cooper 
pairs tunnel due to realistic direction dependent [14] 
|7k P | 2 . This triplet part in the charge Joseph- 
son tunneling process may effectively provide spin 
Josephson tunneling. The rate of change of spin 
S z = E p [ c p t c pT - c pi c pl] due t0 tunneling, i.e., 
S z — i[Ht 1 S z ], gives rise to spin- Josephson tunnel- 
ing current Ij(t) = — 2 3[exp(— 2ieVt)^ ct (eV)] where 



3 



& s (iuJ rs 



Jq dr e iuJmT (T t B(t)A(0)) with B(t) 



E kp T kp [4t( t ) c pt( t ) - c ia( T ) c Pi( T )]- We thus find 
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where (rjk x Vp)j represents j-th spatial component of 

m x t) p . 

We assume that the non-centrosymmetric axes are par- 
allel to the interface of the junction as shown in Fig. 1. 
The probability of tunneling [14] will be the most along 
the direction transverse to the interface. In terms of nor- 
mal state electronic tunneling conductance Gn — In/V 
with //v being normal state tunneling current, we find dc 
(V = 0) charge and spin Josephson currents: 

I j = {jpf^j I sin ^ {.9+Ai + g-A 2 ) + costp ( g+ A 3 + g-M)] 

(8) 

Ij = {j^j 6k sin £ [ cos ^ +^2)- sin ^(^3 + ^4)] , 

(9) 

where ip = 6l — Or, g± = <5cos£ ± 1, 6 £ (0,1) 
is a parameter depending on the model [15] of tun- 
neling matrix element, i.e., the angular dependence of 
Tk p , cos£ = fiL ■ fiR where both ul and tir arc in 
the plane of the interface, and k is the projection of 
n l x fiR along perpendicular to the interface. Further 
— it 5Zaa' v x v \' (1 + ^AA'cos^) describes the depen- 
dence of normal state tunneling current on 5, Ai = 
r ++ + cos(6» L - 6 R ), Ai = T_+ cos 6 L + T+_ cos 9 R , 

A3 = T sin(^L— Or), and A4 =T h sin6*L— T H sin^^ 

with 
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at zero temperature, where K is the elliptic function of 
first kind. 

Charge (8) and spin (9) dc Josephson currents depend 
on both sine and cosine of the global phase difference 
tp between two superconductors. Instead of continuous 
relative phase difference between superconductors in each 
spin-split band, there is a possibility of phase locking 
such that $l, 9r = nir (n G Z). In this case, A3 = 
Ai = and hence Ij(V — 0) oc sin-0 and Ij(V — 0) oc 
cos?/'. This situation is also true for 6l = Or, even in the 
absence of their locking at mr. The critical charge and 
spin Josephson currents are modulated with the angle £. 

For general values of 0l and Or, charge (8) and spin 
(9) Josephson currents may be parametrized as Ij — 
Jf sm(if) + \i) + J^siniip + X2) an d Ij = J s cos(V' + 
Xi), where J^eT/G N ) = CScost J%(eT/G N ) = V, 
J s {e 2 T/G N ) = C8kwo.£, X i = tan"^^^) and X2 - 
tan-^^JE^) with C = [(A, + A 2 ) 2 + (A 3 + A) 2 ] 1/2 , 



and V = [(At - A 2 f + (A3 - Aa) 2 ] 1 / 2 . With the ap- 
plication of external flux, current and phase relationship 
[2, 3] for Ij may be found out. This will determine Jf 
and J| and thereby 8 and £. These determinations will 
predict the value of /}. The spin current may be tested 
if the proposed set-up acts as a source of spin current. 

The spin- Josephson current (9) is proportional to sin £ 
which means is maximum when the orientation of the 
axes of non-centrosymmetry between left and right side of 
the superconductors are transverse, and it vanishes when 
the axes are parallel. This particular property should 
be useful to control spin-Josephson current by orient- 
ing the axes about which the inversion symmetries are 
lost. Therefore it will be useful to build a Josephson spin 
switch. 

To summarize, we have shown that the Josephson 
tunneling process between two superconductors without 
inversion symmetry consists of tunneling of both spin- 
singlet and spin-triplet Cooper pairs. Because of the di- 
rection dependent tunneling matrix element in real junc- 
tion and the momentum dependent phases of the triplet 
components of the gap function, viz. Ak,|t and Ak,xi> 
are opposite, there is a net spin-Josephson current. This 
spin current is proportional to sine of the angle £ between 
the axes of non-centrosymmetry in the bulk superconduc- 
tors in each side of the junction. It may be used to make 
a Josephson spin switch with the controlling parameter 
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